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ABSTRACT 

We consider a lattice model of a semiflexible homopolymer chain in a bad solvent. 
Beside the temperature T, this model is described by (i) a curvature energy Sh, represent- 
ing the stiffness of the chain (ii) a nearest-neighbour attractive energy £„, representing the 
solvent (iii) the monomer density p — where and fl denote respectively the number 
of monomers and the number of lattice sites. This model is a simplified view of the protein 
folding problem, which encompasses the geometrical competition between secondary struc- 
tures (the curvature term modelling helix formation) and the global compactness (modeled 
here by the attractive energy), but contains no side chain information. By allowing the 
monomer density p to depart from unity one has made a first (albeit naive) step to include 
the role of the water. In previous analytical studies we considered only the (fully compact) 
case p = 1, and found a first order freezing transition towards a crystalline ground state 
(also called the native state in the protein literature). 

In this paper, we extend this calculation to the description of both compact and non- 
compact phases. The analysis is done first at a mean-field level. We then find that the 
transition from the high temperature swollen coil state to the crystalline ground state is a 
two-step process for which: (i) there is first a 9 collapse transition towards a compact "liq- 
uid" globule, and (ii) at low temperature, this "liquid" globule undergoes a discontinuous 
freezing transition. The mean-field value of the 9 collapse temperature is found to be inde- 
pendent of the curvature energy Sh- This mean-field analysis is improved by a variational 
bound, which confirms the independence of the 9 collapse temperature with respect to 8^. 
This result is confirmed by a Monte Carlo simulation, although with a much lower value of 
the 9 temperature. This lowering of the collapse transition allows the possibility (for large 
Sfi) of a direct first order freezing transition, from a swollen coil to the crystalline ground 
state. For small values of Sh, the mean-field two-step mechanism remains valid. 

In the protein folding problem, the "liquid" compact phase is likely to be related to 
the "molten globule" phase. The properties of this model system thus suggest that, even 
though side chain disordering is not taken into account, disordering of the backbone of a 
protein may still be a sufficient mechanism to drive the system from the native state into 
the molten globule state. 
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1. INTRODUCTION 

The thermodynamic equihbrium states of small globular proteins in solution depend 
on solvent conditions, including: temperature, pH, salt concentration and denaturant con- 
centration. 

Roughly speaking, there are two ways to go from a fully unfolded swollen coil state to 
the native state: through a single first order folding transition, or through an intermediate 
collapsed phase [1]. The nature of the transition from an unfolded state to this phase, 
which may be generically called a molten globule [2,3], depends in detail on the specific 
protein and solvent conditions, but often has the character of a continuous ^-collapse. The 
subsequent transition to the native state is generally first order in character. Note that, 
in a finite system such as a protein, the use of the expression "phase transition" should be 
interpreted as "strongly cooperative phenomenon." 

In order to elucidate the balance of enthalpic and entropic terms leading to this phase 
diagram we have previously represented the protein folding process using a highly simplified 
homopolymer model, originally due to Flory [4] . In this model the tendency to form local 
order which may have 1-or 2-dimensional character (a helices or f3 sheets) is complemented 
by a tendency to global compactness representing the effect of hydrophobic forces. The 
properties of this model were previously studied numerically by Kolinsky et al [5]. We 
have previously given an analytical mean field treatment of the model [6,7] which was 
restricted to the fully compact phases and yielded a first order freezing transition towards 
a native-like structure. 

In this paper we extend this model by (i) allowing for the presence of vacancies (ii) 
assigning a finite attractive energy between (nonconsecutive) neighbouring monomers. 
This generalization enables us to study the full phase diagram of the model by the method 
of mean field theory, including the high temperature disordered coil state. From a more 
(bio)chemical point of view, points (i) and (ii) above are a first step towards a proper 
treatment of the solvent properties. 

We report elsewhere the results of numerical studies which yield information about 
the structural character of the partially folded states of this model [8]. 

Besides its application for proteins, this model is clearly of interest to describe the 
melting or crystallization of semi-fiexible polymers [4] : in this case, intermediate compact 
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phases are likely to possess orientational order [9,10]. 

The organization of this paper is as follows. In section 2, we recall for completeness 
our previous study of fully compact phases, through the formalism of Hamiltonian paths. 
The method is applied to the case of 1-dimensional structures which may be thought of as 
helix-like (sheets may be treated in a similar way). The extended model, also restricted 
to the case of helix-like structures, is presented at a mean-field level in section 3, together 
with its resulting phase diagram. In section 4, we derive an upper bound for the free 
energy, which strongly suggests that the collapse temperature is indeed independent of 
the curvature energy. Using previous determinations of the 6 temperature, we obtain a 
modified phase diagram, which is then checked numerically by Monte Carlo simulations. 
Connections with other approaches are considered in the conclusion. 

2. SECONDARY STRUCTURES IN A COMPACT PHASE 

In this section, we shall describe how one can study, in terms of a very simplified and 
schematic homopolymer model the thermodynamics of secondary structures in a compact 
phase. Before doing so, let us reconsider the problem of the collapsed chain on a lattice, 
initially without introducing any tendency to secondary structure formation (i.e. for zero 
stiffness). 

2.1 Hamiltonian paths and the compact phase 

We consider a (i-dimensional hypercubic lattice of N = L'^ sites, and its associated 
Hamiltonian paths representing the conformations of a polymer with — 1 links. We 
recall that a Hamiltonian path (HP) is a path which visits all sites of the lattice once 
and only once . In this approach, to compute the entropy of the collapsed but disordered 
phase means to evaluate the number of Hamiltonian paths on the lattice. For simplicity, 
we consider closed paths, but, as is well known in polymer theory [11], (i.e. in the large N 
limit), boundary conditions play only a subdominant role in terms of the free energy. 

Thus we define the partition function Z by: 




(2.1) 



{HP} 



Note that Z is nothing but the total number of HP. One may rewrite (2.1) as: 




(2.2.a) 
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where 

= \Y.^{^ (A.>-0"V(r (2.2.6) 

rr' 

In eq.(2.2), (p{r) is a n— component field attached to each point r of the lattice, and the 
operator Afj? / is 1 if r and r' are nearest-neighbour, and otherwise. The denominator 
is a simple Gaussian integral, equal to exp(— ^ logdet Afsp/). It is equal to one in the limit 
n — > 0, and will be omitted henceforth. 

To prove the identity between (2.2) and (2.1), namely that (2.2) indeed counts all HP, 
we use Wick's theorem. At each site r, one attaches a factor ^(p ^ (r), which can be viewed 
graphically as one line coming into , and one line coming out of site r. The elementary 
contraction (see eq. (2.2.b)) reads: 



if) ■ (fif3 (r ') = SaisAfp / (2.3) 

where a and /? denote any of the n— components of the field 0{r). Since eq.(2.3) conserves 
the index a, each connected part of the path yields a factor n through the summation 
over a. To deal with the one chain problem (or with a single connected path), one must 
therefore extract the term proportionnal to n, hence the factor lim^^o ~ in cq.(2.2.a). By 
doing so [9], paths with more than one connected part, which have higher powers of n, are 
cancelled, which proves the identity of (2.1) and (2.2). 

The exact evaluation of (2.2) is out of reach, and we restrict ourselves to a simple 
approximation, namely the Saddle-Point Method (SPM). Minimizing the exponent of the 
numerator of (2. 2. a) with respect to 0{r) yields: 

^, v(n 

If we assume periodic boundary conditions on the lattice, it is legitimate to look for a 
homogeneous solution to (2.4). We further break the 0{n) symmetry by choosing <^ in a 
given direction, say 1, and a take the solution as: 

(^(r) = (vp,0,...) (2.5.a) 

where <^ is a constant, independent of space. Using (2.4) we get (p^ = 2q, where q = 2d is 
the coordination number of the lattice. In the SPM, the partition function (2.4) reads: 
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where e = exp(l) = 2.718 This result [12] is in excellent agreement with numerical 

simulations and exact results, on various lattices, even in low dimensions. This is due to 
the fact that fluctuations around the SPM are indeed small [2,12]. 

Equation (2.5.b) shows that in the limit of zero stiffness the compact phase has a 
sizeable conformational entropy S/N = log(g/e) , and thus no definite three-dimensional 
shape: a polymer globule (compact phase) is a statistical mixture of an exponentially large 
number of conformations. (We remind the reader that these results are strongly dependent 
on the boundary conditions, assumed here to be periodic). To sum up, we have modeled 
the collapsed phase of a single polymer chain below the 9 point by lattice Hamiltonian 
paths. In this approach, one looses track of the relation between the spatial degrees of 
freedom and the curvilinear abscissa along the chain. 

2.2) Secondary structure in the compact phase 

Following Flory [4] , we may think of each link of the Hamiltonian path as representing 
a helical turn. Since hydrogen-bonds have a tendency to favor long helices, that is to align 
the links of our model, we attribute an energy penalty Eh to the breaking of an helix, that is 
whenever the Hamiltonian path makes a turn (corner). This model of a homopolymer with 
finite stiffness has attracted a lot of attention in the theory of polymer melting [6,9,10]. 
The partition function of the system, at inverse temperature P = j^, reads 



where {HP} denotes the ensemble of all Hamiltonian paths, and NciTi.} denotes the 
number of corners present in path Ti. Following the previous section, we introduce on each 
lattice site r and for each direction a = 1, d an n-component real field 0^. (^) • The 
partition function Z can be rewritten: 




(2.6) 



{HP} 




Z = lim — 

n — >-0 n 




(2.7.a) 



with 




(2.7.6) 
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where the operator A~ , is 1 if r and r ' are nearest-neighbour in direction a, and 
otherwise. Note again that we have omitted the denominator in (2.7.a), since it goes to 
unity when n ^ 0. 

In order to prove the equivalence of (2.6) and (2.7), we will use again Wick's theorem: 
we define, from (2.7.b), the elementary contraction: 



^i"^ (r) • (f ') = 5uv5a^^%- ' (2.8) 

where u and v are component indices of 0ct if) , running from 1 to n. 

Expanding the product over {f} in (2. 7. a), we must choose at each site f either a term 
of the form ^(^^ (r) or one of the form e~^^^0a (0 ' 0-y {'^ ^ corresponding respectively to 
a path going straight through r in direction a, or to a path making a turn at r, from 
direction a to direction 7. By contracting the fields according to (2.8), we construct a sum 
over all self avoiding compact closed paths with appropriate weights, with an additional 
factor n (due to the summation over component index u) for each closed loop. As above, 
we extract the single chain contribution by taking the limit n — 0. This concludes the 
proof. 

The mean field equations (saddle point on (2.7.a)) read: 

7i ^Ea'^a'W + ^e Ey^a' ^a' (r) • (r) 

Assuming periodic boundary conditions, the isotropic homogenous mean- field solution 
takes the form (p(^{f) — (p for any a — 1, d and r. From (2.9), we obtain: 

^^'^ = \A (2.10.a) 
At this mean-field level, the free energy per site reads: 

f=-T\og(^) (2.10.6) 



with 

q{(3) ^2 + 2{d-l) e-^^^ (2.10.c) 
and e = 2.71828... . Note that q{P) plays the role of an effective coordination number . 
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The phase transition of the model described by eq.(2.7) has been studied in detail in 
reference [6] . The ground state consists of a bundle-of-sticks like conformation taking the 
form of straight paths which make turns on the surface. This solution, which has / = 
in the large limit, is not obtained from the saddle point equation (2.9), which is not 
surprising for a structure which is essentially one dimensional. It must nevertheless be 
considered, since it is clearly the lowest energy state of the chain. Within our approxima- 
tions, a first order phase transition must then occur at a temperature Tp for which / = 0, 
yielding an effective coordination number q{P) equal to e. For d = S, Tp = 0.58 Eh- Below 
this temperature, the free energy remains equal to zero. 

Physically, there is a competition between the entropy gain of making turns, and the 
corresponding energy loss. At high temperature, the corners are mobile in the bulk, leading 
to a liquid like structure, whereas at low temperature, the system is frozen, in stretched 
walks, with the corners expelled on the surface. The two regions are separated by a phase 
transition at Tp. 

The average length of a helix is given by 

where U {(3) = — ^log Z is the internal energy. Just above the freezing point, in d = 3, the 
average secondary structure persistence length is equal to ip = 3.78, and is of 0{£ = N^/^) 
in the low temperature phase. Note that in this very simplified picture £f corresponds to 
a typical number of amino-acids per a-helix of the order of 15, since one link corresponds 
to a helical turn, that is 3.6 amino acids. 

Since the density of corners is extensive (i.e. 0{N) for T > Tp), and drops to zero 
(i.e. C(iVi-i/'^)) for T < Tp, the freezing transition is a first order in the thermodynamic 
limit. 

This model can be easily extended to incorporate the existence of /^-sheets [7]: the 
physics of the latter model is very similar to that of the previous model (2.2), namely, a 
liquid-like high temperature phase with no definite /3-sheet structure, and a low tempera- 
ture frozen phase, consisting of stacks of parallel /3-sheets. 
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3. VACANCIES AND GEOMETRICAL FRUSTRATION: A MEAN-FIELD 
TREATMENT 

We wish now to explore the full phase diagram of the protein chain, including its 
disordered (denaturated) coil state. For the sake of simplicity, wc will generalize the 
approach of section 2 to the case of helix-like structures (again, the extension to the (5- 
sheet case is straightforward). We consider a d-dimensional hypercubic lattice of O = L*^ 
sites, and a polymer chain of N monomers, each monomer being at a lattice site. Following 
the preceding discussion, we attribute an energy penalty for each corner (i.e., when two 
consecutive links are not aligned). Furthermore, we now introduce a finite Van der Waals 
attractive energy Sv for each pair of non consecutive monomers separated by one lattice 
spacing. We emphasize once more that this model aims at a simplified description of the 
solvent conditions. 

The partition function reads 



e-^^''^={^^>e-^^^-e^^^''5]."'^-^^-'"--' 5(^n^--iV) (3.1) 

(n^=0,l) {Un} r 

where Uf- denotes the occupation number of lattice site r, and {Hj^} denotes now a self 
avoiding path of monomers on the lattice. The constant factor q-^P^^ in cq.(3.1) ensures 
that one does not count the attractive energy between consecutive monomers. Following 
section 2, it is easy to reformulate this problem by introducing at each site f of the lattice 
(i) n-components fields 0ot (r)' with a = l,....,d (ii) a scalar field if^ff needed to perform 
the Hubbard-Stratanovich transformation on the occupation numbers Uf. Introducing a 
fugacity K per monomer, the grand canonical partition function Z{K) may be written: 



-IP d, 

Z{K) = Ika^- Yl dV'r n ^'^^ ^'^"^ n + Ke-^^-eV^^^-S^.j (3.2.a) 

r a=l r 



with 



and 



r,r 



a.—\ f,r ' 



Bp 



a=l 



l<a<7<d 



(3.2.6) 



(3.2.c) 



it is easy to check, using the approach of section 2, that Zjss in (3.1) is the coefficient of 
in eq.(3.2); we may alternatively write 



91ogZ(K) 



(3.3) 



Our strategy is now straightforward: we evaluate (3.2) via an isotropic saddle point method 
and periodic boundary conditions {(pa.{f) = (p\ '^p — '(/')• With these hypothesis, eq.(3.2.a) 
gives 



logz(i^) = n 



-^(^V'' + \f) + log(l + Ke-^^^'eV^^S) 



(3.4) 



where B is now space independent and given by (3.2.c). Breaking the 0{n) symmetry as 
in section 2, the saddle point equations reads 



i) _ (d^^) K./^ e-^^-eV^'^ (1 + (f^ - 1) e-^^'') 



{2d) 



D 



(1 + (J-1) c-'^^") 



D 



with 



D = l + 



Ke-^^^ev^^ (1 + (d - 1) e"^^'^) 



Equation (3.3) in turn gives 



N = Q 



{dcp"^) Ke-f^'-eV^^ (1 + (t^ - 1) e"'^^'') 



D 



Introducing the monomer density p = we have from equations (3.5) 



(3.5.a) 



(3.5.6) 



(3.5.c) 



{3.5.d) 



^ =dP 



= 2d-\/ps^ p 



leading to a free energy per monomer 



f ^1 ^ fZ{K)\ 
/=--logZ^ =--log^-^j 



(3.6.a) 

(3.6.6) 
(3.6.C) 



(3.7.a) 
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We finally obtain 

/ = -Tlog (^^^ + T^^^ log(l - p) + - de,p (3.7.6) 

This equation is reminiscent of a Flory-Huggins formula [13] for a single chain: the 
first term represents the entropy of the collapsed chain (see eq.(2.10.b)), the second term 
accounts for the entropy of vacancies, and the last two terms represent the attractive 
energy. Note that, we work here with the free energy per monomer rather than with the 
free energy per site. 

As discussed in section 2, the fully stretched (or crystalline) phase makes turns at the 
"surface" , among the vacancies; this phase is not obtained in the saddle point approxima- 
tion. Its free energy per monomer is easily obtained as 

g = -{d-l)e^ (3.8) 

since is the energy between two neighbouring non consecutive monomers. Note that 
for density of vacancies 0, i.e. p = 1, the analysis of section 2 is recovered: one may 
also argue that the limit — oo in equation (3.7.b) implies that the phase with p = 1 
corresponds to the minimum of the free energy /. 

Since our study encompasses the study of the high temperature phase as well, it is 
also interesting to consider the value of the free energy / for the limit p = 0, that we call 
/o. We get 

/o = -Tlog?(/3) + £„ (3.9) 

that is the free energy of the swollen coil, which corresponds, in our approximations, to 
the entropy of a Brownian walk on a lattice with temperature dependent connectivity. 

We now consider the phase diagram implied by the study of equations (3.7.b) (3.8), 
and (3.9). 

Broadly speaking, one should compare the free energies per monomer / and g. A first 
order freezing of the type described previously will occur when g < f. We have chosen to 
present our result in the (t, x) plane, with t = — and x = —. Beside the extremal points 
p = and p = 1, we have to consider the solution p = p* with 

'df 



. . =0 (3.10.a) 



11 



yielding 



log (1 - p*) 



+l+/-=0 



(3.10.6) 



( a) the 9 collapse transition 
From eq.(3.10.b), one may get a non-zero density p*, below the x-independent 9 
temperature tg = 2(i, i.e. tg = & iov d = 3. 
( b) the freezing transition 
One has to compare /(p*) and ^, where p* is given by eq.(3.10.b). It is easy to 
see that for x — > oo, we have at the transition p^ = 0.5, and too = (4 iog'2-2) ' yielding 
too — 3.88 for d = 3. For a; — > 0, we recover the result of section 2, namely pg = 1 and 



Note that in the mean-field treatment of our model, it is not possible to obtain a 
direct transition from the high temperature phase to the crystalline ground state. This 
transition would occur for fo = g, and the resulting temperature is lower than to, as seen 
from eqs(3.8) and (3.9). The fact that one has to go through a dense intermediate phase 
can be traced back to the vacancy entropy in eq (3.7.b). The resulting phase diagram is 
shown in Figure 1. We now study corrections to the mean-field picture. 

4. VARIATIONAL BOUNDS AND MONTE CARLO SIMULATIONS: THE 
MODIFIED PHASE DIAGRAM 

The 9 collapse transition is a 2nd order transition, and is thus not very accurately 
described by mean-field methods, whereas the freezing transition being first-order may be 
quite well described within this framework. For example, we find too high a value of the 
X -independent collapse temperature {te = 6 for d = 3). Numerical estimates [19] give 
t^) ~ 3.64 for x = 0, which may imply a modification of the phase diagram of figure 1. 

A variational bound to the free energy may be easily obtained. Indeed, from equation 
(3.2.c) we have (using Schwarz inequality): 



to = x/log 




, i.e. to 






(4.1.a) 
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where q{P) is the effective connectivity of the lattice, as given by (2.10.c). From (4.1. a) we 
get the bound 



where T) is the partition function of a chain with no curvature energy, with collapse 

energy e^. We first note on this equation that the r.h.s. displays a collapse transition at 
the 9 temperature of the "pure" system, that is the system with no stiffness. This result 
is consistent with the result of section 3. 

As mentioned above, the "pure" 9 point has been calculated numerically in three 
dimensions, leading to a temperature much lower than the mean-field estimate (to ~ 3.64 
vs. a mean-field estimate of = 6). This effect is partly due to the fact that the 9 point 
is a second order phase transition, which is dominated by fluctuations, which are not well 
taken into account at the mean-field level. 

On the other hand, the freezing transition is a first order transition, and thus is not 
driven by large fluctuations. The mean-fleld approximation is then much more reliable. 

If we take these two facts seriously, we obtain the phase diagram represented on flgure 
2. One qualitative modiflcation of the phase diagram is the possibility of a direct freezing 
transition, without passing through a 9 collapse. Indeed, deflning the critical value Xc as 
the intersection between the line of the 9 collapse {t — 3.64) and the freezing curve, we see 
that for: 

X < Xc- the transition to the native (crystalline) state proceeds through a 9 collapse, followed 



by a freezing transition. There is a whole region in temperature, where the system is 
collapsed, but with no deflnite secondary structure, and this may be identifled with a 
"molten globule" state. 



X > Xc- the freezing transition occurs at a higher temperature than the 9 collapse, and thus 
the system will directly freeze, from the swollen coil to the native state. 
Using equations (3.7.b), (3.8), and (3.10) together with the three dimensional numerical 
value of te, we flnd Xc — 15. 

We have checked these predictions numerically by performing Monte Carlo simula- 
tions. We have used a Monte Carlo growth method [20] , which has proven both fast and 
reliable for calculations of thermodynamical properties of random chains. We grow pop- 




(4.1.6) 
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ulations of M = 20000 chains of length N = 200, for various values oi t = T/sy and 
X = Sh/^v There is of course a degradation of the statistics in the frozen phase, and we 
have thus restricted the parameter a; to 3. The numerical results are plotted on figure 2. 
Each calculated point is represented by a star. We see that the agreement is quite good, 
and seems to confirm the theoretical phase diagram predicted above. 

5. CONCLUSIONS 

We have studied a simplified homopolymer model to describe the interplay between 
compactness and local order in a protein chain. In a mean-field treatment of this model, 
the native state can be reached only after a 9 collapse transition towards a liquid globule 
(the collapse transition temperature being independent of the stiffness (or curvature) en- 
ergy). In this globule, helix-like secondary structures progressively grow and ultimately 
freeze. These results are not capable of representing the traditional two-state coil <-> native 
renaturation transition for globular proteins. 

We have argued that this disagreement may be partially explained by the mean-field 
treatment. Using more reliable estimates of the three dimensional tg temperature [19], we 
have used variational bounds, together with Monte Carlo calculations, to show that (i) 
tg is indeed independent of (ii) the phase diagram may be qualitatively modified for 
large enough values of x = Sh/sy, allowing for a unique transition from the coil towards 
the native state. In our model, we get a very high value for the critical Xc, separating the 
two transitions regime from the unique "coil- native state renaturation transition" (xc — 15). 
This value, which shows that a large portion of the phase diagram is occupied by a compact 
liquid or molten globule phase, is certainly not realistic for globular proteins. 

Finally, we think that dealing with an homopolymer model may be inappropriate to 
describe the renaturation transition of globular proteins. The latter has been modeled in a 
simplified theoretical approach by Finkelstein and Shakhnovich [14,15] and by Alonso, Dill 
and Stigter [16]. Transitions of this type have also been studied analytically in a random 
heteropolymer model [17] . Such two-state transitions have their origin in the heteropolymer 
nature of proteins for which hydrophobic residues become segregated within the core of 
the protein, while hydrophilic residues are in contact with solvent. 

Our results show that the homopolymer model has a molten globule like thermody- 
namic phase. Although the model has no side chains and is hence incapable of including 
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the effects of side chain entropy on the native <-> molten globule transition, it does reveal 
another entropic contribution which has not been considered in this context by previous 
workers: the entropy of backbone disordering. 

Thus, even though side chain rotational entropy is undoubtedly important in the 
molten globule state, we have shown that backbone disordering may be capable of driving 
a first order transition, and hence that side chain entropy need not necessarily be the 
driving force for the transition. Thus backbone disordering provides a complementary 
mechanism to the discussion of Bromberg and Dill [18], who argue that the onset of side 
chain disordering may not be sufficiently sharp to drive the native <-> molten globule 
transition. 

Last but not least, it is clear that our results pertain to the case of semi- flexible 
polymers. The intermediate compact phase is then likely to posess orientational (e.g. 
nematic) order [9,10]. 
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FIGURE CAPTIONS 

Figure 1 : Mean field phase diagram in the {t, x) plane. The solid line denotes the contin- 
uous 9 transition from the coil (C) to the liquid globule (G). The dashed line denotes the 
first order freezing transition to the native state (N). 

Figure 2 : Numerical three dimensional phase diagram in the (t, x) plane. The solid line 
denotes the continuous 9 transition from the coil (C) to the liquid globule (G). The dashed 
line denotes the first order freezing transition to the native state (N). The stars show the 
results of the Monte Carlo calculations. The multicritical point corresponds to {tc ^ 3.6, 
Xc — 15). 
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